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ABSTRACT 

This paper investigates GRB 050802, one of the best examples of a Swift gamma- 
ray burst afterglow that shows a break in the X-ray lightcurve, while the optical 
counterpart decays as a single power-law. This burst has an optically bright afterglow 
of 16.5 magnitude, detected throughout the 170- 650 nm spectral range of the UVOT 
on-board Swift. Observations began with the XRT and UVOT telescopes 286 s after 
the initial trigger and continued for 1.2 x 10 6 s. The X-ray lightcurve consists of three 
power-law segments: a rise until 420 s, followed by a slow decay with — 0.63 ± 0.03 
until 5000 s, after which, the lightcurve decays faster with a slope of = 1.59 ±0.03. 
The optical lightcurve decays as a single power-law with qq =0.82 ±0.03 throughout 
the observation. The X-ray data on their own are consistent with the break at 5000 s 
being due to the end of energy injection. Modelling the optical to X-ray spectral 
energy distribution, we find that the optical afterglow can not be produced by the 
same component as the X-ray emission at late times, ruling out a single component 
afterglow. We therefore considered two-component jet models and find that the X-ray 
and optical emission is best reproduced by a model in which both components are 
energy injected for the duration of the observed afterglow and the X-ray break at 
5000 s is due to a jet break in the narrow component. This bright, well-observed burst 
is likely a guide for interpreting the surprising finding of Swift that bursts seldom 
display achromatic jet breaks. 
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1 INTRODUCTION 

Gamma Ray Bursts (GRBs) are the most energetic explo- 
sions that take place in our Universe, with a typical energy of 
10 51 — 10 53 ergs released on a timescale of between a millisec- 
ond and a few thousand seconds. The release of such a con- 
siderable amount of energy ov er such a short period req uires 
an outflow that is relativistic (Mcszaros & Rccs 1997a) and 
is likely to be anisotropic (|Sari et all Il999h . The energy 
within the outflow is released primarily though shocks. The 
GRB is thought to be produced through intern al shocks from 
inter actions between successive shells of ejecta (|Sari fc Piranl 
1 19971 ). Subsequently an afterglow is emitted as the outflow 
is decelerate d through collisionless shocks with the exter- 
nal medium l|Sari|[l997l ). A forward shock propagates into 
the external medium and emits from X-ray to radio wave- 



lengths, and a reverse shock travels back through the ejecta 
and peaks at long er wavelengths than the forward shock 
l|Zhanget alj|2005h . 

The afterglow can reveal many properties of the progen- 
itor an d its surroundings . Currently, the most up to date 
model l|Zhang et alj [2005) allows the use of the temporal 
and spectral indices to indicate the nature of the surround- 
ing medium, whether it is a uniform density medium, or 
a medium with a density that is decreasing radially as ex- 
pected for a stellar wind e.g n tx r~ 2 , where n is the particle 
density and r is the radius from the GRB. The indices also 
indicate the location of the observing band relative to the 
synchrotron self absorption frequency v a , the peak frequency 
v m and cooling frequency v c . 

Some GRB afterglows have a period of sl ow decline in 
their lightcurves (F(t) tx t~ a , where a ~ 0.5) (|Zhang et al.l 
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l2005h . which is generally accepted to b e due to continued 
energy injection (|Rees fc Mesz aros 1998). This e nergy could 
be du e to a central engine that is long lasting ([Dai fc Lu.l 
Il998al ; IZhang fc Meszaros.1 1200 ll ). later shells catching up 
and colliding with slow er shells that were emitted earlier 
(R ees fc Meszaros! [r998), or the slow release of energy store d 
in the form of Poynting Flux l|Zhang fc Kobavashil I2005T ). 
The end of energy injection is signalled by an increase in 
decay rates of the afterglow lightcurves to F(t) oc t~ a , where 
a ~ 1. 

Currently the most effective observatory for studying 
GRBs and their early afterglows is Swift, which has now 
been in operation for over two years. It has the ability to 
observe emission ranging from 7-rays to optical with the 
three o n-board telescope s, namely the Burst Alert Telescope 
(BAT; iBarthelmv et ail 2005), the X-ray Telescope (XRT; 
iBurrows et alj 2005) and the Ultra Violet and Optical Tele- 
scope (UVOT; [Roming et all 2005). The BAT covers a large 
area of the sky allowing of order 100 bursts to be detected 
per year. When the BAT has been triggered by the GRB, 
the satellite slews automatically allowing the XRT and the 
UVOT to begin observing as soon as possible, usually within 
100 s of the burst trigger. Swift is thus able to observe both 
the initial 7-ray explosion and the following early afterglow. 

This paper looks at GRB 050802, a burst that appears 
to challenge the standard picture. At the end of the shallow 
decay the X-ray afterglow breaks to a steeper decay, while 
the optical afterglow continues to decay as a power-law with- 
out a break. If both the X-ray and optical arise from the 
same component and if the X-ray break were due to the end 
of energy injection then the optical lightcurve would be ex- 
pected to break at the same time. Such be haviour has been 
noted recently in 6 Swift GRB afterglows l|Panaitescu et all 
2006a) of which GRB 050802 is currently the best example. 
This burst has a X-ray lightcurve with two distinct breaks 
and it was observed for 1.2 x 10 6 s. It also had an optically 
bright afterglow (~16 magnitude at early times) that was 
well sampled in 6 filters of the UVOT up to 1 x 10 s s and 
thereafter observed with the white UVOT filter. In this pa- 
per we discuss the possible models that could explain this be- 
haviour. We will use the convention flux F — t~ a with 
a and /3 being the temporal and spectral indices respectively. 
We assume the Hubble parameter Ho = 70kms _1 Mpc -1 
and density parameters Qa = 0.7 and Q m =0.3. Uncertain- 
ties are quoted at la unless otherwise stated. 



2 SWIFT AND GROUND BASED 
OBSERVATIONS 

The BAT was trigger ed by GRB 0508 02 at 10:08:02 UT on 
the 2 nd August 2005 (jBand et alj|2005h ■ The lightcurve rises 
for 5 s to the first of three peaks and has a T90 = (13 ± 2) s 
(90 % confidence level). The fluence in the 15 - 350 keV 
band is (2.8 ± 0.1) x 10~ 6 ergcm~ 2 to 90% confidence 
level (|Palmer et al.ll2005l ). Observations with the XRT and 
UVOT began 289 s and 286 s respectively, afte r the BAT 
trigger (|Band et alJliooBI : iMcGowan et aLlfeoOSal ). Both the 
XRT and UVOT continued to observe until 1.2 x 10 6 s after 
the burst trigger. 

The XRT began observations by locating the burst with 
Image Mode (IM). After the burst was located, data were 



taken in Windowed Timing (WT) mode for 163 s. A fad- 
ing unca talogued source w as found within 8 of the BAT 
position |Band et al.ll2005i ) and was confirmed as the X-ray 
counterpart of GRB 050802. 480 s after the burst trigger, the 
XRT changed modes and continued observations in Photon 
Counting (PC) mode. 

UVOT observations showed a fading, uncatalogued 
source at RA = U h 37 m 05.69 s , Dec = 27° 47 12.2" 
l|McGowan et al.l (2005a). Following the trigger, a se- 
ries of automated exposures were taken in the three optical 
and three ultra-violet filters. A log of the observations is 
given in Table [T] The observations consisted of an initial 
100 s 'finding chart' exposure in the V-band, 10 s exposures 
in each passband for 7 rotations of the filter wheel, followed 
by a sequence of ~ 100 s and ~ 900 s exposures. Later 
observations (after 1 x 10 s), were taken in the UVOT 
white filter. 

The afterglow was also imaged and detected with 
the 2.6 m Shain Telescope, 8hrs after the burst, with 
R and i' band magn itudes of 20.6 and 20.2 respectively 
(|Pavlenko et al.ll2005h . Spectroscopic observations were car- 
ried out with ALFOSC on the Nordic Optical Telescope. 
Several absorption features were de tected provid ing a red- 
shift m easurement of z = 1.71 |Fvnbo et al.l 2005a & 
2005b). iFvnbo et all (|2005aD discovered an extended source 
within 1" of the UVOT afterglow location, which they pro- 
pose could be the host galaxy. 



3 XRT AND UVOT DATA REDUCTION. 
3.1 X-ray Data 

The XRT data were reduced using the XRT pipeline soft- 
ware V2.0. Source and background counts were taken from 
the cleaned event files using extraction regions, in order to 
construct spectra and lightcurves. Events of grade 0-2 were 
used for the WT mode and 0-12 were used for the PC mode. 
For the WT mode data we used a 40 pixel strip for the source 
extraction region and a 40 pixel strip for the background ex- 
traction region. In the PC mode data the first 2.5 ks of data 
were found to be piled-up and so required the use of an annu- 
lar extraction region. The size of the region that was affected 
by pile up was determined by comparing the radial profile of 
the afterglow with a model of the XRT point spread function 
(PSF). The radial profile and the model PSF are inconsis- 
tent at radii less than 2.5 pixels (6"), so we used an inner 
radius of 2.5 pixels and an outer radius of 30 pixels (71"). 
The size of the source extraction region for the rest of the PC 
data was 30 pixel (71") radius. For all PC data, background 
counts were extracted from a circular region of radius 80 pix- 
els (189"). Appropriate response matrices (RMs) were taken 
from the Swift calibration database, CALDB 20060424 and 
effective area files were constructed using the standard XRT 
software. A correction factor was calculated and applied to 
the piled-up section of the XRT lightcurve to account for 
the excluded, piled up pixels. Fortunately, the source was 
not located near the bad pixel columns and so no correction 
was required. 
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3.2 UVOT Data 

The UVOT event files were screened for bad times (e.g. 
South Atlantic Anomaly passage, Earth-limb avoidance) and 
the images were corrected for Mod-8 noise using the stan- 
dard UVOT software. The images were transformed to sky 
coordinates and then corrected for the ~ 5" uncertainty in 
the aspect of the spacecraft pointing using bespoke software. 

Counts for the afterglow were extracted using an aper- 
ture of radius 4" for the optical filters and 5" for the ultra- 
violet filters. Background subtraction was performed using 
counts extracted from a larger region offset from the source 
position. The measured count rates were aperture-corrected 
to radii of 6" for the V, B and U filters and 12" for the UV 
filters. These were then translated to magnitudes using the 
standard UVOT zero-points (see Table [lj. 

The optical lightcurve is shown in Fig.[T] The lightcurve 
from each filter of the UVOT was normalized to the V filter. 
The normalization factor for each filter (see Table [2} was 
calculated by taking the average count rate through the filter 
in the 400s - 1000s time range, and dividing this value by 
the average count rate in the V filter over the same time 
interval. Later observations were obtained with the UVOT 
white filter; these were normalized to the equivalent V count 
rates as follows. The optical-UV spectral energy distribution 
(SED) was modelled using the average count rates from the 
V, B, U, UVW1, UVM2 and UVW2 filters in the 400s- 
1000s time range. The optical/UV response matrices were 
then used to predict the ratio of the V to white count rates. 



3.3 Combined X-ray and UV/optical Spectral 
Energy Distributions 

SEDs were produced spanning the optical to X-ray range for 
early (400s - 1000s) and late (35ks - 55ks) times. For each pe- 
riod, the average count rates of the exposures in each UVOT 
filter were used to produce the optical spectral values. For 
the X-ray part of each SED a spectrum was extracted in the 
relevant time range. 



4 RESULTS 

4.1 X-ray and Optical Lightcurves 

The 0.2- lOkeV X-ray lightcurve of GRB 050802 is shown in 
Fig. [2] A visual inspection shows a complex behaviour with 
an initial rise followed by a flat and then a more rapid de- 
cay. The X-ray lightcurve was first modelled using a broken 
power-law. The best fit parameters were ot\ = 0.55 ± 0.03, 
a 2 = 1.59 ± 0.03 and break time 4600 ± 260 s. However the 
X 2 /D.O.F = 81/57 corresponds to a null hypothesis proba- 
bility of only 0.02, and the model systematically deviates 
from the observed lightcurve at the earliest times. Hence, 
the lightcurve was modelled using a double broken power- 
law (i.e. a model with 3 power-law segments). This provides 
a better fit with x 2 /D.O.F = 64/55; according to the F-test 
the 3-segment power-law fit gives an improvement at the 3<r 
confidence level with respect to the 2-segment model. The 
values of the best fit parameters are shown in Table [3] 

In the best fit model, the X-ray lightcurve first rises 
with a slope ai = -O.SO+^g until 420 ± 40 s. At this 



point, the lightcurve breaks for the first time and a shal- 
low decay begins with a 2 = 0.63 ± 0.03. This phase ends at 
5000 ± 300 s when the lightcurve starts to decay steeply with 
q 3 = 1-59 ±0.03. 

To look for spectral variations over the course of the 
decay, we split the lightcurve into soft (0.2-2keV) and 
hard (2 - 10 keV) X-ray lightcurves. The soft and hard X- 
ray lightcurves are shown in Fig. [3] where we also show the 
softness ratio, which we define as (Cs — Ch)/(Cs + Ch) 
where Cs is the soft X-ray count rate and Ch is the hard 
X-ray count rate. A constant can be fitted well through the 
softness ratio time series (x 2 /D.O.F = 30/29) revealing that 
there is no significant spectral evolution with time. 

The optical lightcurve is shown in Fig. [T] and it is well 
fitted with a single power-law decay with a temporal index 
Qi =0.82 ±0.03 (x 2 / D -O.F = 71/63). The lightcurve has no 
obvious colour evolution within the wavelength range of the 
UVOT over the duration of the burst afterglow. 

iFvnbo et all (2005a, 2005b) observed an extended 
source within an arcsecond of the afterglow location, 1.5 days 
after the burst trigger. They suggested that it might be the 
host galaxy and they provided a combined magnitude for the 
afterglow and extended source of R = 22.5. We have looked 
for the host galaxy by coadding the UVOT white filter im- 
ages taken at late times after the afterglow has faded be- 
yond detection. The UVOT white filter observations were 
taken until 1.2 x 10 6 s and are shown in Fig. [1] The sum- 
mation of the white exposures provides a deep 3cr upper limit 
of 23.4 mag, which is equivalent to V = 23.5 mag for the af- 
terglow spectrum. From these observations, we can confirm 
that there is no significant contribution from the host galaxy 
to the lightcurve in the UVOT spectral range for the first 
60 ks, while the afterglow is still detected by the UVOT. 

To determine the earliest time at which the optical 
lightcurve could have broken, a 2-segment power-law was 
fitted with the second segment decay rate set to be the same 
as the third segment in the X-ray lightcurve. We determined 
the 3<7 lower limit for the break time by adjusting the time 
of the break until we obtained Ax 2 = 9 with respect to the 
single power-law fit. The lower limit to the break time was 
found to be 19 ks after the burst trigger, significantly later 
than the second X-ray break. 

In Fig. [4] we show the X-ray/optical ratio, which we 
define as (Cx — Co)/(Cx + Co) where Co is the opti- 
cal/UV count rate normalized to the V filter and Cx is 
the X-ray count rate. Initially, we tried fitting a con- 
stant across the entire time range. This provided a poor 
fit with x 2 /D.O.F = 168/62, implying that the changes in 
the X-ray/optical ratio are highly significant. To investi- 
gate the behaviour of the X-ray/optical ratio we fit a func- 
tion to it before and after the second X-ray break time of 
5000 s. The first 5000 s was fit by a linear relationship be- 
tween (Cx — C'o)/{Cx + Co) and log(t) with a gradient of 
0.24 ± 0.01 (x 2 /D.O.F = 47/47). After 5000 s, the best fit 
linear relationship has a x 2 /D.O.F = 7/12 and gradient of 
-0.09 ±0.01. The evolving X-ray/optical ratio implies that 
the optical/UV to X-ray SED is changing throughout the 
afterglow. 
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Table 1. UVOT observations of GRB 050802 given in (aperture corrected) count rates and magnitudes. Count rates were aperture 
corrected to 6" for V, B and U filters and to 12" for UVWl, UVM2, UVW2 niters to enable the use of the UVOT zero points to convert 
the count rates to magnitudes. The zero points are V = 17.83X0.09, B = 19.12X0.12, U = 18.34X0.23, UVWl = 17.82X0.02 ,UVM2 
= 17.19 ±0.23, UVW2 = 17.82X0.02, White = 19.78X0.02. 



4.2 X-ray and Optical Spectra 

The results of the X-ray spectral analysis are shown in 
Table [3] All spectra were fitted using an absorbed power- 
law model. The absorption component includes both pho- 
toelectric absorption from our Galaxy and from the host 
galaxy of the GRB. The Galactic column density was fixed 



at Nh = 1.78 x 10 cm 2 (jDickev fc LockmarJll99un and 
the column density at z = 1.71 was allowed to vary. The 
spectral slopes and column densities measured for the 3 seg- 
ments of the X-ray afterglow show no evidence of evolution. 
There is evidence for absorption from the host galaxy of 
the GRB in each segment; the column density is consistent 
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Filter 


Ratio 


B 


0.39 


U 


0.45 


UVW1 


1.71 


UVM2 


14.46 


UVW2 


12.45 


White 


0.13 



Table 2. For 6 UVOT filters, the normalization factor required to convert to V count rate. For B, U, UVW1, UVM2 and UVW2 filters, 
the normalization factor was calculated by taking the average count rate through the filter in the 400s - 1000s range, and dividing this 
value by the average count rate in the V filter over the same time interval. The average count rates were then used to create an optical/UV 
SED, which was used to predict the ratio of White to V count rate. 
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Figure 1. Optical Lightcurve of GRB 050802 in 4 of the 6 filters available. The UVM2 and UVW2 filters were excluded from all plots as 
the data was not constraining. The count rates inside the 10 s exposures, between 400 s and 1000 s were summed and averaged for each 
filter and the rest of the exposures were normalized against this value. The normalized count rate in each exposure was then divided by 
the average value for the V filter. Upper limits are given to 3cr. 



between the 3 segments of the X-ray lightcurve, with an 
average value of 2.9 x 10 21 cm~ 2 . 

Optical/UV to X-ray SEDs in the time intervals 400s- 
1000s and 35ks - 55ks were created as described in Section 
3.3. For each SED, a power-law fit accounting for Galac- 
tic and GRB host-galaxy dust and photoelectric absorption 
was applied (see Fig. [5}. For extinction in our Galaxy, a 
fixed dust component was used with E(B- V) = 0.03, using 
the Milky Way extinction curve (He3[l992). For extinction 
in the host galaxy the fit was tried using the SMC extinc- 
tion curve and then with the Milky Way extinction curve. 
The fitting was applied to the 400s - 1000s SED as these data 
have better signal to noise than at late times. The fit using 
the SMC extinction curve produced j3 = 0.79 ± 0.02 and 



X 2 /D.O.F = 134/104. The fit using the Milky Way extinction 
curve returned (3 = 0.86 ± 0.02 and xV D -O.F = 120/104. 
We also tested a model in which a cooling break resides in 
between the optical and X-ray bands. In this case, a broken 
power-law model was tested in which the spectral indices 
have a fixed difference of A/3 = 0.5. For consistency, the 
model was tried using the SMC and MW extinction curves. 
The broken power law using the SMC extinction curve re- 
turned E break = 0.010lg;gg|keV and j3 a = 0.89 ± 0.01 with 
X 2 /D.O.F = 125/103. The fit using the Milky Way extinc- 
tion curve gave E break = 0.004jl o ;g u | and = 0.89 ± 0.04 
with x 2 /D.O.F = 119/103. Overall, the fits with the Milky 
Way extinction curve provide the best x 2 7 1 D.O.F and for 
this extinction curve there is no significant improvement to 
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Segment 


Time at which slope breaks(s) 


Energy index (3 


Ar H (10 20 cm- 2 ) 


Temporal Index a 


1 (initial rise) 

2 (shallow decay) 

3 (steep decay) 


420 ± 40 
5000 ± 300 


0.87 ± 0.08 
0.89 ± 0.04 
0.88 ± 0.04 


28 ± 10 
31 ± 5 
28 ± 5 


-o.801° ;S 

0.63± 0.03 
1.59 ± 0.03 



Table 3. Spectral and temporal analysis of the X-ray lightcurve of GRB 050802 fitted with 3 segments. The value of Njj refers to the 
absorption at z = 1.71. 




Time (s) 

Figure 2. X-ray Lightcurve of GRB 050802. A double broken power-law is found to fit the X-ray lightcurve with breaks at 420 ± 40 s 
and at 5000 ± 310 s. 



the fit by replacing a power-law with a broken power-law. 
The model parameters from the different model fits and the 
implied total (Galactic and GRB host galaxy) extinction in 
the UVOT bands are given in Table [4] 

Unusually, we find that the Milky Way extinction 
curve best fits the SED of GRB 050802. In compari son, 
no other burs t s in t he samples of ISchadv et al.l (|200?t ) or 
IStarling et~al1 (|2006h . which consist of 7 Swift GRBs and 10 
BeppoSAX GRBs respectively, are fitted best with a Milky 
Way extinction curve. Since the extinction curve is unusual, 
we determined the Nh /Ay to see how this compares to other 
GRBs. For GRB 050802, the N H /A V ratio was found to 
be 4.5 (±2.3) x 10 21 . The mean G RB N H /A V ratio fo r a 
MW extinction law in the sample of ISchadv et al. | (|2007h is 
i.7tii x 1 21 , so GRB 050802 is consistent with the mean 
Nh/A v in ISchadv etafl (|2007h to within la. This implies 
that the ratio of dust and gas surrounding GRB 050802, is 
fairly typical for GRBs. 

Because we have shown that the absorption does not 
change significantly with time, we do not expect the extinc- 
tion to change either. Therefore, in fitting the 35ks - 55ks 



SED we froze the extinction and absorption at the best fit 
values found for the 400s- 1000s SED. This resulted in a fit 
with (3 = 0.99 ±0.02 (x 2 /B\G- F = 27/15), which implies 
a null hypothesis of only 3%. Furthermore, the value of f3 
obtained in this fit is inconsistent at 99 per cent confidence 
with j3z = 0.88± 0.04, the spectral index of the third seg- 
ment of the X-ray lightcurve. If we repeat the fit to the 
35ks-55ks SED with a fixed spectral index of /3 = 0.88, we 
obtain a \ 2 /D.O.F = 35/16, which implies the model is 
rejected at 99 per cent confidence. 



5 DISCUSSION 

In summary, we have determined that the X-ray lightcurve 
can be divided into 3 segments with two breaks at 420 s 
and 5000 s. Throughout the afterglow, there is no evidence 
for X-ray spectral evolution. We also find that the opti- 
cal lightcurve decays as a single power-law, which has a 
slightly steeper decay than the second segment of the X-ray 
lightcurve. We have determined the earliest possible time 
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Figure 3. Top panel: The soft (0.2-2keV) and hard (2-10keV) X-ray lightcurves. Bottom panel: The softness ratio, defined as (C's — 
Ch)/(Cs +Ch), where Cr is the hard X-ray count rate and Cg is the soft X-ray count rate. There is no evidence for spectral evolution 
within this time period as can be seen by the constant ratio. 



Model Parameters — Models for 400s - 1000s SED - 35ks - 55ks SED 



& Host Extinction 


Power-law 


Power-law 


Broken power-law 


Broken power-law 


Power-law 




MW extinction 


SMC extinction 


MW extinction 


SMC extinction 


MW extinction 





0.86 ±0.02 


0.79 ±0.02 


0.89 ±0.04 


0.89 ± 0.01 


0.99 ±0.02 


Energy break (keV) 






004+ ' 005 
u.uu^_ 003 


010+ 009 
u.uiu_ 008 




E(B - V) host 


0.18 ±0.02 


0.09 ±0.01 


0.18 ±0.02 


0.10 ±0.02 


0.18* 


A v (Mag) 


26 ±4 


20±4 


29 ±6 


30 ±4 


26* 


1.6±0.2 


0.8±0.1 


1.6±0.2 


0.8±0.1 


1.6* 


A B (Mag) 


1.4±0.2 


1.2±0.1 


1.5±0.2 


1.1±0.1 


1.4* 


A v (Mag) 


1.7±0.2 


1.4±0.2 


1.8±0.2 


1.4±0.2 


1.7* 


Auvwi (Mag) 


2.5±0.3 


1.8±0.2 


2.6±0.3 


1.8±0.3 


2.5* 


AuvM2 (Mag) 


2.8±0.3 


1.9±0.2 


2.9±0.3 


1.9±0.3 


2.8* 


Auvw2 (Mag) 


2.8±0.3 


1.9±0.2 


2.9±0.3 


1.9±0.3 


2.8* 


X 2 /D-O.F 


120/104 


134/104 


119/103 


125/103 


27/15 



Table 4. The model parameters were determined from the fitting of a power-law and separately a broken power-law to the 400s- 1000s 
SED; the fits were repeated for the MW and SMC extinction curves. Also shown, is the power-law with MW extinction fit to the 35ks- 
55ks SED. Parameters marked with * arc fixed at the best fit values found for the power-law fit to the 400s -1000 s using MW extinction. 
All models use a fixed value of the Galactic extinction E(B — V) = 0.03 and Galactic absorption Nu = 1.78 X 10 20 cm" 2 . The observed 
host extinction in each filter is provided for each fit; note that these were derived from the fit parameter E(B — V)host an d were not fit 
parameters in their own right. 
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Figure 4. Top panel: the X-ray and optical lightcurves. The optical light curve is shown for 4 out of the 6 filters. The count rates in 
each filter have been normalized to the V filter, as discussed the caption to Fig. [T] The X-ray lightcurve has been binned according to 
the length of each exposure in the optical lightcurve. Different symbols refer to the filters used for the optical observation in each time 
bin. Bottom panel: the ratio between the X-ray (Cx) an d optical (Go) lightcurves. 



that the optical lightcurve could have changed its decay rate 
to match that of the third X-ray segment to be 19 ks, which 
is 14 ks after the last X-ray break. 

In the following subsections, we will look in detail at the 
X-ray afterglow to examine the origin of the X-ray breaks 
and to determine the best fitting closure relations between 
the temporal and spectral slopes. Then we shall discuss the 
question of why the X-ray lightcurve breaks at 5000 s while 
the optical lightcurve continues to decay as a power-law. 
Several mechanisms will then be investigated to find a com- 
patible model for the production of the X-ray and optical 
afterglow. 



5.1 The X-ray Afterglow 

We begin by looking at the first break to see what we 
may learn. The break occurs at 420 ± 40 s, which is pre- 
ceded by a rise from 286 s with a slope of -0.80±°;^. The 
break is then followed by a gradual decay with a slope of 
a — 0.63 ± 0.03. The rise and the successive slow decay 
can be regarded as a broad peak in the lightcurve and it 
may be attributed to 'flaring' activity (jBurrows et al . 2005b; 



iNousek et al.l [20 06) or could be the early phase of the jet 
interacting with the external mediu m, giving rise to the for- 
ward shock (|Kumar fc P iran 2000aj). However, it is not pos- 
sible to discriminate between the two mechanisms. Flares 
often have different spectra to the afterglow, but if this is a 
flare it may be so small we would not expect to see any spec- 
tral variability. On the other hand, we can not tell if this is 
the early rise of the afterglow as we have too little data prior 
to the peak. An additional contribution to the second seg- 
ment from a large flare may also be excluded as a power-law 
fit to the X-ray lightcurve between 1000s and 3000s results 
in a slope of a = 0.64 ± 0.05 with x 2 / D -O.F = 21/38. 

Typically, X-ray afterglo ws are analysed by applying 
the closure relations given bv lZhang et alj (|2005) to the X- 
ray spectral and temporal indices. These relations are an 
aid in determining the location of the observed X-ray band 
relative to the synchrotron frequencies (y a , v m , u c ) and the 
environment in which the the burst occurs. 

We skip the first segment of the X-ray lightcurve as 
the closure relations are not applicable to a lightcurve that 
is rising. We begin by applying the simplest non-injected 
relations to the second segment, in which «2 = 0.63 ± 0.03 
and /?2 = 0.89 ± 0.04, the only closure relation that agrees 
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Figure 5. Optical and X-ray SEDs between 400s and 1000s (top) and 35ks- 55ks (bottom). The top panel shows the fit using SMC dust 
extinction (grey; offset by a factor of 1.1) and Milky Way extinction (black). The first 6 points in each panel show the 6 filters of the 
UVOT ranging from V to UVW2. Points after 0.2 keV show the X-ray spectrum from 0.2-10keV. In the top and bottom panels, the 
dahsed line represents the weighted mean of the spectral indices found in the 3 segments of the X-ray lightcurve, /3 = 0.88. 



without the requirement of energy injection at 3<r confidence 
level, is 



3/3 



(1) 



This relation applies to a wind or ISM medium with vx > 
v c , but, this relation is not appropriate for GRB 050802 as 
the value of the energy index p, determined through j3 = p/2 
is p= 1.78 ±0.08. Values of p < 2 require a diffe r ent se t of 
closure relations provided by IZhang &T Meszarosl (|2004 ) , of 



which none satisfies the second decay of GRB 050802. More 
complex closure relations, which include a thi rd parameter 
to acc ount for energy injection, are provided bv lZhang et al.l 
(2005). In these models the luminosity is assumed to evolve 
as L(t) = Lo(t/tb)~ q where q is the luminosity index and tb is 
the time for the formation of a self-similar solution, which is 
approximately the deceleration time. The luminosity index 
q must be < 1 for the energy injection to affect the afterglow 
dynamics. We found that there are two injected closure re- 



10 



(2) 



lations that are consistent with 02 and fa. The first relation 
is 

q , (2 + <?)/? 
Q= 2 + ^2~ 
This closure relation is for a wind medium with slow cooling 
electrons with either vx < v m or v m < vx < v c and where 
g= —0.27 ± 0.04. The second possible energy injected rela- 
tion is 

(2 + q)P 



a = (q-l) + 



(3) 



This relation is for a uniform medium with slow cooling 
electrons with either v m < vx or v m < v x < v c and where 
g = 0.51 ± 0.04. For the closure relations given in Equations 
[land El /3 = (p - l)/2 and sop = 2.79 ± 0.08. 

When the closure relations are applied with the indices 
of the third segment, where «3 = 1.59 ± 0.03 and /?3 = 
0.88 ± 0.04, we find that the most consistent closure relation 
without energy injection is 



a = (3/3 + l)/2 



(4) 



which is consistent at 3.5a confidence. All other non-injected 
closure relations are ruled out at a minimum of 4a. Equa- 
tion!?] is for electrons that are slow cooling within the range 
v m < vx < v c in a wind medium. Again a value for the en- 
ergy index may be determined though /? = (p — l)/2 with 
p — 2.76 ±0.08. Generally, the energy index is expected to 
remain constant throughout the X-ray afterglow , unless it is 
changed by energy injection l|Zhang et al.ll2005T ). 

Therefore, the application of the closure relations to the 
X-ray temporal and spectral indices suggests that the sec- 
ond segment is energy injected and the third segment is not, 
and therefore that the break between the second and third 
segments at ~ 5000 s is due to the discontinuation of energy 
injection. At this point it is worth pointing out that the sec- 
ond and third segments of GRB 050802 are consistent with 
the second a nd third segm e nts of the canonical lightcurve 
presented by iNousek et all |2006). Furthermore, from the 
application of the closure relation s above we wo u ld com e to 
the same ph y sical explanation as INousek et al.l (|2006l ) and 
IZhang et al.1 (|2005h for this temporal behaviour: that the 
break between the two segments corresponds to the end of 
energy injection. Having looked at the X-ray afterglow in 
isolation, we will now investigate whether this physical pic- 
ture is compatible with the observed optical emission. The 
standard afterglo w is based on the assump tion of a single rel- 
ativistic outflow l|Meszaros fc Reesll 1997a! ) and we will begin 
by considering models of this type before considering multi- 
component outflows. 

5.2 Can the X-ray and optical afterglow be 
explained by a single jet? 

We start by considering an afterglow in which the emission 
in both the X-ray and optical bands is produced predomi- 
nantly by the forward shock. If we compare the decay of the 
optical lightcurve with the segments of the X-ray lightcurve 
we find that the optical decay is similar, although not iden- 
tical to, the slow 'energy injected' decay of the second seg- 
ment. This would suggest that the optical lightcurve is also 
characteristic of an energy-injected afterglow. Since the X- 
ray lightcurve breaks at 5000 s, suggesting the end of energy 



injection, and the optical lightcurve continues with a shal- 
low decay indicative of energy injection, we first consider the 
possibility that the electrons responsible for producing the 
optical emission continue to receive energy while the X-ray 
emitting electrons do not. 

There are two possible mechanisms for continued en- 
ergy injection. These are the arrival of additional shells of 
material at the shock region or the release of energy stored 
in the form of Poynting flux. In the first case, the injected 
energy changes the fireball dynamics, but can only change 
the balance of X-ray to optical emission and produce a break 
in one of the lightcurves if v m or v c pass through one of these 
bands. However, we can rule out the passage of one of these 
breaks through the X-ray or optical afterglow because there 
is no break in the optical lightcurve and there is no change 
in the X-ray spectrum. This therefore rules out continued 
energy in the form of blast waves and so we consider contin- 
ued energy injection in the form of Poynting flux. Energy is 
supplied in the form of Poynting flux when there is a rotat- 
ing compact object that has a magnetic field, B. The flux is 
converted to kinetic energy when kink instabiliti es cause the 
magnetic field lines to reconn ect in the outflow (|Drenkhahnl 
l2002j ; lGiannios fc Sprurd2006D . The energy density of Poynt- 
ing flux is oc B 2 , and within the reconnection region the en- 
ergy available p er particle for particle acceleration is oc B 
(Giannios 2006). To decrease the photon energy so that fur- 
ther emission is below the X-ray band, requires a decrease 
of the average particle energy. The XRT is sensitive over 
more than a decade in photon energy, therefore a decrease 
in the average photon energy from above the XRT energy 
band to below the XRT energy band requires a factor >10 
change in B and therefore a fact or of > 100 cha nge in the 
energy density of Poynting flux l|Giannios! [20061 ). This im- 
plies a large change in the overall rate of energy injection, 
which is inconsistent with the continued power-law decay in 
the optical. 

A reverse shock is an alternative means for producing 
optical emission without X-ray emission and so we now ex- 
amine the possibility that this could sustain the optical emis- 
sion of GRB 050802 after the X-ray lightcurve has broken. 
The reverse shock will travel back through the ejected mat- 
ter and will emit as long as it is passing through ejecta. 
The reverse shock e mission will cease when it passes over 
the last ejected shell (|Zhang fc Meszaroslliool '). The shock 
crossing time is expected to be approximately the duration 
of the GRB, but if energy continues to be injected into the 
afterglow by the arrival of further shells the reverse shock 
may continue for as long as the energy injection phase. If 
the end of energy injection occurs at 5000 s, as suggested by 
the X-ray lightcurve, then we should expect that the reverse 
shock must cease at approximately 5000 s. Therefore, this is 
not the required mechanism as the optical lightcurve has no 
break until at least 19 ks. 

iPanaitescu et all (|2006aT ) have proposed that it is pos- 
sible to sustain a power-law decay in the optical, while the 
X-ray lightcurve shows a break due to the end of energy in- 
jection, if the microphysical parameters are changing with 
time. The microphysical parameters that must change to 
keep the optical temporal decay as a power-law are: the 
fraction of post-shock energy in the magnetic fields es, the 
fraction of post-shock energy that is given to electrons ej 
and the blast-wave kinetic energy E (No usek et al.l 120061 : 
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iPanaitescu et al1l2006al ). These parameters are expected to 
evolve with Loren tz factor V, if they evolve at all. In the sce- 
nario proposed bv lPanaitescu et al.l (|2006al ) . the cooling fre- 
quency u c lies between the optical and X-ray bands. In this 
case, there should be a slope change of A/3 = 0.5 between 
the X-ray and optical bands and v c mov es to lower frequen- 
cies with time (|Panaitescu et alj|2006ah . For GRB 050802, 
through the fitting of a broken power-law with a fixed 
A/3 = 0.5 to the 400s- 1000s SED (see Section 4.2), the 
maximum energy, at which a break could be present, was 
determined to be 0.02 keV. 

We examined the optical to X-ray SEDs to see if the X- 
ray and optical lightcurves can be explained by v c moving 
to lower frequenc i es with time as expected in the model of 
IPanaitescu et a l. (2006a). In this model the largest change 
in the ratio (Cx — Co)/{Cx + Co), where Cx is the X-ray 
count rate and Co is the V band count rate corresponds 
to the motion of v c from its highest allowed value to below 
the optical pass band, at which point the optical emission 
lies on the same power-law spectral segment as the X-ray 
emission. In the broken power-law fit to the 400s - 1000s 
SED, the highest allowed value of v c at 3<r confidence is 
0.06 keV. Therefore we refited the 400s- 1000s SED with v c 
fixed at 0.06 keV, then reduced the value of v c to below the 
UVOT spectral range and so determined that the small- 
est {Cx — Co)/(Cx + Co) value allowed by this model is 
-0.10. The observed range of (Cx - Co)/(Cx + Co) ex- 
tends to much lower values than this (Fig. [4]). We find that 
the weighted mean of the data points more than 35 ks after 
the BAT trigger is (Cx - Co)/{Cx + Co) = -0.38 ± 0.08, 
which is > 3<r below t he lowest allowed value of —0.10, thus 
allowing the model of IPanaitescu et al.l (|2006al ) to be ruled 
out. 

We now look in to the optical to X-ray SED further and 
ask: to what extent can single component models be ruled 
out completely? As discussed in Section 4.2, a power-law fit 
to the late time (35ks-55ks) SED provides a poor fit and 
the best fit value of /3 is softer than the X-ray spectral index 
for segment 3, which includes the 35ks-55ks time interval. 
In this respect, we note that there is no evidence for any 
evolution in the X-ray spectral index at any point in the af- 
terglow and the weighted mean of the indices found in the 3 
segments of the X-ray lightcurve gives /3 = 0.88 ± 0.03. To see 
why the fit to the late time SED is poor we show the X-ray 
spectrum extrapolated to lower energies using this spectral 
slope (/3 = 0.88) as a dashed line in both panels of Fig [5] 
The dashed line lies above the optical and UV data points 
as expected for an afterglow with significant extinction. In 
Section 4.2, we determined that at 3<r confidence, the host 
galaxy extinction is greater than 0.8 and 1.1 in the B and 
U bands respectively, no matter which models we choose for 
the extinction and continuum. In addition, the Galactic ex- 
tinction is 0.1 for both the B and the U bands. After correc- 
tion for extinction, the B and U band fluxes in the late time 
SED lie significantly above the dashed line. As described in 
Section 4.2, this power-law can be ruled out with >99 per 
cent confidence, implying at 99 per cent confidence that the 
region responsible for the X-ray emission can not produce 
all of the optical emission at late times, no matter where v c 
and v m lie. Since single component outflows are unable to 
represent the late time SED satisfactorily, we now consider 



multi-component outflows to see which of these models are 
able to reproduce the observed lightcurves and SEDs. 

5.3 Multi-component Outflows 

A multi-component outflow is one consisting of two or more 
components that have different bulk Lorentz factors. The 
simplest model is one of two components: a narrow jet 
and a second wider, but slower jet that surrounds the nar- 
row component. This geometry may be generated by the 
narrow jet giving rise to a wider and slower component 
(|Kumar fc Gra not 2003), where the bulk F decreases over 
time and produces a jet with angular structure. Alterna- 
tively, the two components may be formed at the same 
time for example when neutrons and protons decouple in 
a neutron-rich, hydroydynamically accelerated jet from a 
neutron star or from th e neutron-rich accret io n disk of a 
collap sed massive star (|Vlahakis et all l2005h . IPeng et al.l 
l|2005h show that within such a neutron-rich jet, the wide 
component has greater energy than the narrow component 
so that the wider component is able to dominate at late 
times. In both cases, the narrow component will produce 
the X-ray emission, whereas the optical emission will be pro- 
duced by the wider component as it travels at lower Lorentz 
factors than the narrow component. This picture can be used 
to explain bursts such as GRB 021004 and GRB 030329 
jPeng et al.ll2005h . The wider jet is expected to reach the 
undisturbed medium at a slightly later time than the nar- 
row jet. Thus, the wide component may be observed as a 
later start for the optical lightcurve when compared to the 
X-ray lightcurve, or on an established optical lightcurve, the 
arrival of the wide jet to the external medium may be shown 
as a late rise/bump or a change in the decay slope. 

There are three options within a two-component jet we 
will investigate to see if one could have produced the ob- 
served lightcurves of GRB 050802. All of these options in- 
volve energy injection, which is introduced through shells of 
material. The options are: energy injection with the ejected 
shells wide enough to refresh both components of the jet 
with a Lorentz factor that decreases over time (so that after 
5000 s the shells no longer reach the shock in the narrow 
component), the cessation of energy injection at 5000s in 
one component only and, finally, the continuation of energy 
injection for the duration of the afterglow, but with the nar- 
row component producing a jet break at 5000 s. 

The first option is that energy injection is distributed 
uniformly throughout both jets and continues for a sustained 
period of time. Initially, the shells have a wide enough angu- 
lar distribution, and travel at a F large enough, to reach the 
forward shock produced by both components. Both the X- 
ray and the optical afterglows will appear energy injected. 
As time goes by, the Lorentz factor of the freshly ejected 
shells decreases. Eventually, the ejected shells will not travel 
at a Lorentz factor large enough to reach the shocked region 
in the narrow component, but will still be able to reach the 
wide component. This will cause a break to be seen only 
in the X-ray afterglow. However, the freshly ejected shells 
will still travel at F large enough to reach the wide com- 
ponent and therefore the optical lightcurve will continue 
to appear energy injected. This much appears to be con- 
sistent with the data because, as discussed in Section 5.2, 
the optical lightcurve appears to be energy injected for its 
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entire observed duration and does not break with the X- 
ray lightcurve. However, as the wide component is slower 
than the narrow component, the wide component will absorb 
shells (i.e receive energy) at a faster rate than the narrow 
component. Therefore, the optical emission is expected to 
decay less rapidly than the X-ray emission. For GRB 050802, 
the opposite is actually observed: in the optical the decay 
index ao ~ 0.82 ± 0.03, while the X-ray second segment 
has a decay index of «2 = 0.63 ± 0.03, thus this model is 
excluded. 

The second option calls for the X-ray break at 5000 s 
to have been caused by the cessation of energy injection in 
the narrow component, while injection continues in the wide 
component. However, this implies that the energy is injected 
as a hollow jet at late times and we consider this to be an 
unattractive solution on physical grounds, and lacking in 
testable predictions. 

The third option is that energy is continuously in- 
jected into both components and the X-ray break at 
5000 s results from a jet b reak of the narrow component 
l|Panaitescu fc Kumadl2004h . The spectral indices of the X- 
ray segments support this model as there is no spectral 
change over the break, as is expected for a jet break. A 
jet break may be produced within a laterally expanding or 
a non-expandin g jet. The closure relatio ns for these cases 
are provided by IPanaitescu et alj (|2006bl ). In order for this 
scenario to be viable we have to find a post jet-break clo- 
sure relation that has consistent energy injection, external 
medium, synchrotron frequencies and electron spectral in- 
dex to that of the pre jet-break closure rel ation. The closure 
relation that satisfies this uses Eq 35 of IPanaitescu et all 
|2006al ), and is: 

a =(l + 2/?)-|(l- g )(/? + 2) (5) 

This relation is independent of the GRB host medium and 
is for a laterally expanding jet with v m <v< v a . The rela- 
tion provides a value of q = 0.38 ± 0.04, which is consistent at 
the 99 per cent confidence level with the value of q obtained 
through the closure relation in Equation [3] The consistency 
with Equation implies that this burst originated in a uni- 
form ISM and not a wind medium as suggested in Section 
5.1. 

Furthermore, we can use the second X-ray break time 
to infer the opening angle of the narrow beam and to put 
a lower limit on the angle of the wide component. We 
calculated the a ngles using the following expression from 
iFrail et all (|200ll ): 

e ict = 0.057 t/ et j [^ r ) (htJ (6) 

where Ei ao is the isotropic energy of the burst, n is the 
density of the medium and r\ is the efficiency of convert- 
ing energy in the ejecta into 7-rays. We use values of 
E iso = 0.2 x 10 53 ergs (in the range 15-350 keV) and GRB ef- 
ficiency 77 = 0.2. Approximating the circumburst medium as 
a spherical, uniform gas cloud of radius R, the particle den- 
sity n, may be estimated simply as n = Nb / R- The cloud ra- 
dius for gas an d dust to survive t he GR B, may be estimated 
using Fig. 8 of |P crna fc Lazzatil (|2002h . This shows that the 
obscuring medium must extend to R ~ 10 20 cm, implying 
a number density of n ~30 cm -3 . Thus, the value for the 



narrow component is estimated to be 6j e t,narrow ~ 1° and 
for the wide component the opening angle is #j e t,wide ^ 8°. 
We recognise that these angles are only approximate be- 
cause they depend weakly on the bolometric energy, the 
efficiency of the burst and on the particle density of the 
external medium. 

This model requires that energy be continually injected 
into both components of the jet for the observed duration 
of the burst. The time scale over which this energy injection 
takes place is more than an order of magnitude larger than 
10 3 — 10 4 s, which is the duration of energ y injection required 
to explain the canonical X-ray afterglow llNousek et all 2006 ; 
I Zhang et al.l [20051; ICapalbi et allkoOfj ; iTashiro et aTThood ; 
iKrimm et all 120061 ). This situation may seem unusual, but 
GRB 050802 is not the first burst that has required dis- 
crete or continuous energy injection for an extended period. 
A few examples of bursts that have required long dura- 
tion energy in jection are GRB 021004, XRF 050406 and 
GRB 060729 (IBjornsson et all |2004 iRomano et"all 120061 : 
iGrupe et al . 2006). The continuation of energy injection well 
beyond the jet break time will produce shallower post jet 
break slopes than the non-injected case and these may be 
mistaken for the discontinuation of energy injection. This 
would lead to the surprisi ngly small number of r eported jet 
breaks in the Swift era ( Willing ale et all [20061 . However, 
the mechanism that could be responsible for long duration 
(> 10 4 s) energy injection is uncertain. Currently, the most 
favoured scenario is one in which the duration of the central 
engine is short, but where the ejecta are em itted with shells 
of varying velocities (|Rees fc Meszaroslll998l ). The shells will 
order themselves before reaching the external medium and 
will form a continuous flow of shells. In this case, small time 
scale fluctuations need not be observed in the power-law 
lightcurves and the central engine is not required to be par- 
ticularly long lived. 

This model would allow the break to occur only in the 
X-ray lightcurve because the narrow and wide components 
travel at different F, and the narrow component is beamed 
into a smaller angle. Therefore, the two component outflow 
with continued energy injection and a jet break in the nar- 
row component is our favoured hypothesis for producing the 
observed lightcurve of GRB 050802. 

The result from this paper has wider implications for 
our understanding of GRBs. Segments two and three of the 
X-ray lightcurve of GRB 050802 have the same character- 
istics as the segment s two and three of t he canonical X-ray 
afterglow lightcurve l|Nousek et alj|2006f ). If we had looked 
at the X-ray data alone, then our explanation for the second 
and third X-ray segments of GRB 050802 wo uld have been 
equiva lent to the explanation proposed by iNousek et all 
( 2006) and lZhang etafl (|2005T ) for the canonical lightcurve: 
an energy injected decay followed by a 'normal' decay. How- 
ever, when the optical observations are taken into account, 
this explanation is no longer attractive. Instead, a more 
plausible explanation is a two component outflow with the 
break in the X-ray lightcurve likely to be due to a jet break. 
Our results suggest that without the optical lightcurve the 
interpretation of GRB afterglows may not be correct. 
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6 CONCLUSION 

This paper has investigated GRB 050802, an unusual burst 
because the optical decays as a power-law, while the X-ray 
breaks twice. We have analysed optical and X-ray data gath- 
ered from the XRT and UVOT instruments onboard Swift. 
Observations continued with both the XRT and the UVOT 
up until 1.2xl0 6 s. 

The analysis of the afterglow began by looking at 
the origin of the three power-law segments in the X- 
ray lightcurve. The X-ray lightcurve commences with a 
rise until 420 s. The second segment decays slowly with 
Q2= 0.63 ±0.03 until 5000 s. The third and last segment 
decays faster with a slope of 03 = 1.59 ± 0.03. The optical 
lightcurve decays as a single power-law with oto = 0.82±0.03. 

Through modelling of the 400s- 1000s X-ray to optical 
SED, it was determined that the best fit is a power-law with 
Milky Way extinction. However, when extinction in the host 
galaxy is accounted for, the optical points lie above the X-ray 
power-law at late times, indicating that the optical afterglow 
could not be produced by the same component as the X-ray 
emission. This result rules out single component afterglows 
with 99 per cent confidence. 

The next logical step was to examine a multi-component 
outflow. Several variations of a two component jet, the sim- 
plest multi-component outflow, were investigated. Of all 
the possibilities, we find the most physically self-consistent 
model is one in which both a narrow and wide component 
are energy injected for the duration of the observed after- 
glow and the X-ray break at 5000 s is due to a jet break in 
the narrow component. 

This paper has wide implications for the GRB commu- 
nity. We have found that the explanation obtained by exam- 
ining only the X-ray data of GRB 050802 differs remarkably 
from the answer obtained by examining the optical and X- 
ray emission together. Our results suggest that without the 
optical lightcurve the correct interpretation of GRB after- 
glows may not be possible. 
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